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Abstract 

We consider Ising model on edge-dual of uncorrelated random networks with arbitrary degree 
distribution. These networks have a finite clustering in the thermodynamic limit. High and low 
temperature expansions of Ising model on the edge-dual of random networks are derived. A detailed 
comparison of the critical behavior of Ising model on scale free random networks and their edge-dual 
is presented. 
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I. INTRODUCTION 



It is evident that a variety of natural and artificial systems can be described in terms 
of complex networks, in which the nodes represent typical units and the edges represent 
interactions between pairs of units jll 0, Qj- Clearly, identifying structural and universal 
'eatures of these networks is the first step in understanding the behavior of these systems 
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Intensive research in recent years has revealed peculiar properties of 
complex networks which were unexpected in the conventional graph theory n|. Among 
these one can refer to the scale free behavior of degree distribution 5], P{k), where degree 
denotes the number of nearest neighbors of a node. From another point of view one 



is interested in the effect of structural properties of 
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behavior of systems living on these networks 
and Ising model (or in general Potts model) are typical examples of statistical mechanics 
which have intensively been studied on uncorrelated random networks with given degree 
distributions Q,E,2i,22,3,24,3,E,22,3,Q- B Y uncorrelated random network we 
mean those in which the degree of two neighbors are independent random variables . These 
networks are identified only by a degree distribution, P{k), and have the maximum possible 
entropy. The locally tree-like nature of these networks provides a good condition to apply 
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the recurrence relations to study the collective behavior of interesting models [22, 
is seen that, depending on the level at which the higher moments of P{k) become infinite, 
one encounters different critical behaviors that could be derived from a landau-Ginzburg 
theory Q]. This in turn reflects the mean field nature of these behaviors. 



In this paper we are going to study the Ising model on the edge-dual of uncorrelated 
random networks with a given degree distribution. These kinds of networks have already 
been introduced in the context of graph theory [l^. Given a network G, its edge-dual G 
can be constructed as follows, see also figure((H): one puts a node in place of each edge of G 
and connects each pair of these new nodes if they are emanating from the same node of G. 
Such networks have been useful among other things for the study of maximum matching 
problem jsij as well as topological phase transitions in random networks 32]. 
The interesting point about these networks is that they have generally a large degree of 
clustering even when the underlying network is tree-like. Due to this high clustering direct 
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FIG. 1: Some part of the edge-dual of a Bethe lattice with z = 3. Empty circles and dotted lines 
denote the nodes and the edges of the edge-dual network respectively. 



study of structural properties of such networks or of physical models denned on them is 
usually very difficult. However one can use this duality to adopt the techniques used in 
the context of random networks (e.g. the generating function formalism 
obtain interesting results for such networks j33j. For example it was shown 
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that the 



edge-dual of a random scale free network with P(k) oc k 7 will be a network whose degree 
distribution behaves like P{k) oc A; -7 for large degrees where 7 = 7 — !. 



Our basic result is depicted in table HI where we have compared the critical behavior of 
Ising model on scale free networks and their edge-dual networks. On the way to this basic 
result we have also studied as a preliminary step the Ising model on the edge-dual of Bethe 
lattices. We have also developed a systematic high and low temperature expansion for the 
Ising model on edge-dual networks for arbitrary degree distributions. 

Moreover as a byproduct we have also shown that there is a simple relation between the 
partition function of an Ising model on a tree in which each spin interacts with its nearest 
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and next nearest neighbors and the partition function of an Ising model on its edge-dual 
with only nearest neighbor interactions but in the presence of magnetic field. 



The paper is organized as follows. In section |H] a relationship between the Ising model on 
a tree and on its edge-dual is derived. High and low temperature expansions of the partition 
function of the Ising model on the edge-dual of a random network are given in section UTTl In 
section HVl we use recurrence method for the study of Ising model on the edge-dual of Bethe 
lattices. The same method is applied in section IS to study the critical behavior of Ising 
model on the edge-dual of scale free networks. The conclusions are presented in section IVI1 

II. A RELATION BETWEEN ISING MODEL ON A TREE AND ITS EDGE- 
DUAL NETWORK 

Let us consider an Ising model (with values of spin taking only ±1) on a tree graph G 
with nearest and next nearest neighbor interactions of strength J\ and J2 in the absence of 
magnetic field. The hamiltonian is 



where < ij >i and < ij >2 denote the nearest and next nearest neighbors respectively. 
For any given configuration of spins we can assign a unique configuration of spin variables 
(again taking values ±1) to the edges of the graph: any edge which connects two nodes 
having the variables Si and Sj is assigned a value S(ij) := SiSj, see figure (J2J). On the other 
hand for any configuration of spins on the edges, S^j), there are two possible configuration 
of spins on the nodes, which are obtained by flipping all the spins Si on the graph. Therefore 
there is a two to one correspondence between the spin configurations on the nodes of the 
graph and the edges of the graph. 

Now if we write the above hamiltonian in terms of spins of edges we get 



where k is the common nearest neighbor of nodes i and j. But this is the hamiltonian 
of an Ising model on G, the edge-dual of the tree, with nearest neighbor interactions of 



E — —J\ ^ SiSj — J2 ^ SiS. 

<ij>i <ij>2 



(1) 



e - - Ji Yl •%) - J 2 S(tk)S(kj), 

<ij>i <ij>2 



(2) 
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strength J\ = J 2 in presence of a magnetic field of magnitude h — J\. 

E^-h^S.-J, £ S-S-, (3) 

i <ij>i 

where i and j now denote the nodes of the edge-dual graph. Taking into account the 
relation between configurations mentioned above we obtain 

Z(J U J 2 , h = 0,T) = 2Z(J X = J 2 , J 2 = 0, h = J u T), (4) 

in which Z and Z are respectively the partition functions of the Ising model on G and 
G and T denotes the temperature. We should stress that the above relation is true only 
for tree graphs, since the presence of loops in G puts constraints on the values of spins 
which are assigned to the nodes of G. That is for any loop in G, the product of spins on its 
edges (equivalently the nodes of G) should be +1. Taking into account all these constraints 
makes the calculation of the partition function very difficult in the general case. 



In the thermodynamic limit the two models have the same free energy density, that is 
f = f, where / := ^ T ^ lZ and / := ~ T ^ Z are respectively the free energy of the Ising model 
on tree and its edge-dual. Here we have set the Boltsmann constant equal to one. Note that 
for a tree graph, N = N — 1. Moreover, it is easy to see that the average magnetization of 
a spin in G 

m = -af, (5) 

is equal to the average correlation of two neighboring spins in G 

< S,S, >= (6) 

Knowing the average magnetization of spins in G as a function of magnetic field, we can 
write the free energy density from equation (JHJ) as follows jsoj 

/( Ji, h, T) = I {m{ti) - l)dh! -,h-h. (7) 

Jh 

Differentiation of the right hand side with respect to h correctly gives the magnetization 
m and the integration constant — J\ — h can be understood from the fact that for very large 
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FIG. 2: To each edge (ij) we can assign a spin which is defined as the product of spins at the end 
point nodes of that edge i.e. := SiSj. 

magnetic fields when the first integral vanishes, the partition function is dominated by the 
configurations where all the spins are up, and hence the free energy per site is given by 
- Ji - h. 

Obviously if we replace J\ with J 2 and h with J\ we obtain /, the free energy per site of 
Ising model on G with nearest and next nearest neighbor interactions. In this way any 
nonanalytic behavior of / will appear in / too. 

III. HIGH AND LOW TEMPERATURE EXPANSIONS OF ISING MODEL ON 
EDGE-DUAL NETWORKS 

Consider a random network which we denote by G. We assume that the probability 
of each node having degree k is equal to P(k). Furthermore we assume that there is no 
correlations between the degrees of adjacent neighbors. The average number of nearest 
neighbors will be denoted by n\ :=< k >= J2kkP(k). Degree distribution of nearest 
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neighbors 2l| is given by n(fc) = ^g 1 . Thus the average number of next nearest neighbors 
will be n 2 :=< k > J2 k (k - l)U(k) =< k 2 > - < k >. 

The edge-dual of such a network is denoted by G. Correspondingly every quantity pertaining 
to the dual network G will be designated by a tilde sign. 

In the following we will study the Ising model with only nearest neighbor interactions on G. 
The partition function is 

Z = £ ^3, ( 8 ) 

where the sum in the exponential is over the nearest neighbors on the edge-dual of G. 
Here we have used the notation K := if. 

A. High Temperature Expansion 

The above partition function can be written in a form appropriate for a high temperature 
expansion j^]. To this end we write the exponential in the form 

e*£<S> S ^ = cosh Z (K) J] (1 + S-S- tanh(K)), (9) 

<i~j> 

where L := ^y 2 - is the number of edges in the G. Inserting this in equation (jHJ) and 
expanding the product we get a series of terms each corresponding to a subgraph of G. 
Summing over spin configuration only terms which represent closed loops will survive 
and we arrive at the following expression for the partition function 



Z = cosh Z (K)2^ V L{C \ (10) 

closed loops c 

where fj = tanh(iT), N = is the number of nodes of G and L(c) is the perimeter (the 
number of edges) of the closed loop c. 

Clearly at high temperatures the first and the second terms corresponding to triangles 
and squares need be kept in the expansion. 

Z = cosh L (K)2*(N A f + N^tf + ....) (11) 
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The number of triangles in G has two parts: first, each triangle of G appears as a triangle 
in G too. Secondly, by definition of the edge-dual network, every triple of edges emanating 
from the same node in G make a triangle in G. The number of the latter types of triangles 
is given by the number of distinct choices of three edges of a node, summed over the nodes 
of G. Thus 



NA = N A + J2( h ), (12) 
i 3 

in which fcj is the degree of node i in G. Since in an uncorrelated random network the 



number of triangles is a finite quantity in the thermodynamic limit [34] , we can neglect the 
first term compared with the second one which has an infinite contribution in this limit. The 
same argument is applicable to the case of squares so we can approximate these numbers by 



#a«£( l ) = ^Y, k (k-l)(k-2)P(k), (13) 

i 3 6 - k 

k- IV 

#o«£( ! ) = ^E^-#-2)(fc-W). 

i 4 4 - k 

Note that these relations become exact in the case of tree structures even for finite N. 



B. Low Temperature Expansion 

We now return to equation (JHJ), the original relation for the partition function. Note that 
we can rewrite it as 

Z= e^ (( ^ 5(lj))2 ~ fcl) , (14) 
{%?)} 

where the first sum in the exponential is over nodes of G and the second is over nearest 
neighbors of that node. Let us write Z in a simpler way 

Z = e- kL £ (15) 
{%?)} 

where L is the number of edges in G and := ^(y)- Using the following identity 
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we find 

Z = e -^(^)^/ 2 ^ J Dxe~^M +2ximi \ (17) 
{%*)} 

where Da; = n« dxi and i runs over all the nodes of G. We can now perform the sum over 
the spin configurations in the integrand. To this end we note in view of the definition of nii 

V e - k ^i Ximi = J2 e~ R ^<^> (Xz+Xj)S{t]) (18) 

{S(ij)} { s (ij)} 

where Y,<ij> sums over all the links of G. The sum can be transformed to 



<ij> 

<ij> 



<ij> 

Putting all these together we find 



, . ) N ' 2 [ Dxe~%^M^>) (1 + e -2K(xi+Xi)y (2Q) 

2n J <ij> 

The product ll<jj>(l + e~ 2K<yXi+x ^) can now be expanded as a series of terms each cor- 
responding to a subgraph g of G. 

For any node i of a the graph G, a factor e ~ 2Kz ^ should be taken into account in which 
Zi is the degree of that node in the subgraph. If a node i dose not belong to the subgraph, 
Zi — 0. Any subgraph determines uniquely the sequence of integers {zf, i — 1, . . . N}. Note 
that Zi < ki Wi. For each such sequence the integral can be easily calculated yielding 

Z = e ~ kL+ ^ ^ k * J2 e~ 2i ^» Zl(k% ~ Zi) • (21) 
g 

It is the central result of this subsection which can be used for a low temperature 
expansion of Ising model on G. This formula incidentally shows that each subgraph g and 
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its complement (the graph obtained when one removes all the links of g from G) give the 
same contribution to the partition function. 

At very low temperatures, K — ► oo, only the empty graph for which all Zi = and its 
complement for which all Zj = ki contribute yielding 

Z = 2e~ kL+ ^ k * 2 = 2 e - kL+ T N E^ 2 P(k) = 2e -iff (*)+f JV(fe 2 > ^2) 
resulting in a free energy per site equal to 

fo = ~((k 2 ) - (k)) (23) 

where we have used the relation K = ^. 
The next to leading order term comes from subgraphs which have only one link, (we multiply 
their contribution by 2 to account for their complements). This will give 

Z = Z {1 + e~ 2ki - k ^- 2) + ■■■) (24) 

<ij> 

where the sum is over all the links of G. This can be written as follows 

Z = Z (l + N^e 4R « e - 2k ^ » + •••) (25) 

where <C3> denotes the average with respect to the two point function P(k,k'), the 
probability of two nodes of degrees k and k' to be neighbors. For uncorrelated networks 
one has P(k,k') = (2 — 5k t k>)n(k)U(k'). This procedure can be followed for higher order 
contributions. 



IV. ISING MODEL ON EDGE-DUAL OF BETHE LATTICES 

The Bethe Lattice j3| is defined cis ct rc gular network where all nodes have the same 
degree z. Let us consider Ising spins on the edges of this network and let them interact with 
an external magnetic field h and with each other if their corresponding edges are incident 
on the same node of the Bethe lattice. We can write the average magnetization of a spin 
lying on an arbitrary edge using the following recurrence relation 
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FIG. 3: The partition function of spins on the branch stems from S can recursively be written in 
terms of the partition functions of branches now stem from nearest neighbors of S. 



; a ( o)-e-y ( o) (26) 

""#) + <-V-(») 

where H := ^ and <?+(0) and <?-(0) are respectively the partition functions for the system 
of spins on one side of the central spin when it is up or down. That is 



{Sr} 



(27) 



where as before K := ^. Note that in this partition function only spins on one side of 
5 (for instance the right hand side spins, denoted by {Sr}) appear. Similar to <?s(0) one 
can define a partition function gs{l) which gives the partition function of the branch of the 
lattice which stems from a node at layer / where the value of its spin has been fixed to 
S. These partition functions can be related to each other recursively as follows, see also 
figure 0: in the right hand of S there are z — 1 spins which interact with S and with each 
other. We can write gs(l) as a sum over different configurations of these spins. For each 
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configuration we will have a term proportional to g r + {l + l)g z T 1 ^ r (l + 1) where r will be the 
number of up spins in such a configuration. Moreover we have to consider another factor 
which takes into account the Boltzmann factor associated to this configuration of spins. 
Energy of a configuration in which r of these spins are up is the sum of three parts, a part 
given by their interaction with the external magnetic field equal to — h(2r — z + 1), a part 
from their interactions with spin S equal to —J(Sr — S(z — 1 — r)) and finally a part given 



by interactions between themselves equal to — J( 
Summing up the above arguments we arrive at 



r(r— 1) 
2 



(z-l-r)(z-l-r-l) 



r(z — 1 — r)). 



2-1 



(/) = V( Z 1 ) e #(2r-, + l) + K(5( 2 r-, + l)-2r(,-l-r)+i^i^i) g r (/ + ^-l-r^ + 1} . (28 ) 



r=0 r 



Returning to equation ([26)1. magnetization of the central spin can be rewritten in a simpler 
form 



e 2H - xl 



e 2H - e~ 2yo 



m 



e 2k + xl e 2k + e- 2 ^ ' 



(29) 



where we have defined xi := =: When the magnetic field is positive we have 

g~{l) < 9+(l) thus yi is a positive quantity which plays a role similar to the magnetic field 
and thus can be interpreted as the local field experienced by a spin at distance I from the 
central spin S. Now using equation (|2Hjl. the recurrence relation for reads 



-In 



( z-i - N 

N e (2i?+2/ !+1 )r+K(«-l-2r-2r(2-l-r)) 

r 



2 — 1 

_ W2fr+ W+ i)r+A-(2r-*+l-2r(2-l-r)) 



E 



V 



(30) 



/ 



Setting H = and starting from distant (/ ^> 1) spins with ?/ < 1 one could obtain 
the values of y for deeper spins in a step by step manner using the above relation until one 
arrives at yo- Equation ()29j) tells us that we will have magnetization in this case only if yo is 
different from zero. It is evident that a stable nonzero solution for yo is possible only when 
the right hand side of the recurrence relation for y t ( when plotted versus yi+i) has a slope 
greater than or equal to 1. The equality will provide the critical temperature of the system 
which turns out to be given by 
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z - 1 



e --^c2r(2-l-r) 2 sinh(X c (2< 



■r - z + 1)) 



r 



1. 



(31) 



z 



1 



Er( 




r 



Unfortunately it is not possible to derive a closed relation for K c . In figure (J2J we have 
computed this quantity numerically and compared it with the corresponding quantity in the 
Bethe lattice itself. In the latter case the critical temperature reads 0| 



As figure (0J) shows T c is much grater than T c which is as expected due to the larger 
number of interactions in the edge-dual network. In the figure we also show a linear fit, 
T c — a,o + a\z, to the numerical data for T c with ao = —3.53 ± 0.02 and a\ = 1.97 ± 0.002. 
As long as the critical behavior of the system is concerned we expect to see a standard mean 
field behavior as in the usual Ising model in spatial dimensions greater than d c = 4. We will 
further discuss these issues in the next section. 

V. ISING MODEL ON EDGE-DUAL OF RANDOM NETWORKS 

In this section we generalize the results of the previous section to the case of Ising model 
on the edge-dual of random networks with a given degree distribution P{k). In this case a 
spin on the edge of such a random network will encounter k\ — 1 and k2 — l nearest neighbors 
at its right and left hand sides respectively. These numbers are random variables given by 
the degree distribution of nearest neighbors in the random network U(k). 

A. General Arguments 

Along the lines of section IIVI we can write the magnetization of a spin on an edge of 
random network with end point nodes of degrees k\ and &2 as 



tanh K c 



1 



(32) 



z 



- 1 



e 2H _ g-yo(fci) -2/0O2) 



(33) 



^fcifc 2 



e 2H _|_ e -yo{k 1 )-y (k 2 ) 
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FIG. 4: T c ( in units of J) for the Ising model on Bethe lattices (lower curve) and their edge- 
dual(upper curve). 

where now y's are random variables depending on distance and degree of end point nodes. 
Obviously magnetization of an arbitrary spin is given by 



m 



n(fc!)n(A; 2 



,2i? 



-yo(ki)-yo(k2) 



ki,k 2 



e 2H _|_ e -yo(ki)-yo(k2) 



(34) 



As before here we used the notation xi{k) := e~ Vl ^ and X\{k) := ^jnx where gs(l', k) is 
the partition function for the cluster beyond the spin S at distance I from the central spin. 
As in the case of Bethe lattices (figure (jHJ)) these quantities can be related to gs'(l + 1; ^')' s 
by the following relations 



9s{l',k) = ^^ e H{2r-k+l)+k{S{2r-k+l)-2r{k~l-r)+ ( ^^^-) £ JJ ^ JJ k b ), 

r=0 r/Ha=l b=r+l 

(35) 

where the second sum is over different selections of r distinct spins from the set of k — 1 
neighboring spins after assigning indices 1 to r to them. Thus the relation for Xi(k) gets the 
form 
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.H2r+K(k-l-2r-2r(k-l-r)) yrk-l-r _ 

or in terms of y's 



W(*0 = -«n(±^ ^k-l-r/k-l - } _ (37) 

p H2r+X(2r-A;+l-2r(A;-l-r)) V , p~Ea=l 

Z^r=0 t—tk— 1— r/fc— 1 ° 

Let us also derive a relation for the average energy of Ising model on the edge-dual of 
random networks in the absence of magnetic field. First note that we can write this quantity 
as a sum over the interaction energy associated to the spins on the edges emanating from 
the same node of G, that is 

£ = E e > = E( E -Js {ij) s m )- (38) 

i i <(ij)(ik)> 

Thus the thermodynamic average of the above quantity reads 

< E>=N"£P(k) < e k >, (39) 

k 

where < e k > is the average energy associated to a node of degree k. We are able to write 
this quantity in terms of g's by summing over different configurations of spins on the edges 
of such a node. As before we sum over configurations in which r spins out of these k spins 
are up. For each such configuration we include an appropriate Boltzmann weight as before. 
Denoting by e k '■— — J( r(r 2 ~ 1 ' > + ( fc ~ r )(^~ r ~ 1 ) _ r (j% _ r )) the interaction energy of these spins 
we obtain 



EJU e k e*^+*^F^-r(*-r)) Er/k nU g+{0 . ka) nUi g _ {0 . kb) 



< > ~ f>rr(r-l) , (fc-r)(fc-r-l) ~. ~ ~ . 7TTT~ ' ^ 4U J 

r/k 



EJU e^^-tM) Er/fc nUl g+{0] ka ) nU <7-(0; 



After some algebra this relation takes the following simpler form in terms of y's 
jf k(k-l) Eto2r(fe-Qe-^ 2r(fc - r) E fc - r / fc e-E::>o(M \ 

Equations (}3*T|) and (jH)) are the exact relations for magnetization, effective fields y 
and energy of the system. 
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B. The Effective Medium Approximation 



In this section we simplify the relations obtained in the previous subsection using the 
effective medium approximation 0, 3 applied satisfactorily to the study of Ising model 
on uncorrelated random networks. It is believed that this approximation takes in a good 
way into account the effects of high degree nodes which play an essential role in determining 
the critical behavior of the system specifically in inhomogeneous network having scale free 
degree distribution. 

To this end we rewrite the relations derived above as if y's are independent of k, the degree 
of the end pint nodes. This is achieved if we use the same g for all the spins which are at 
the same distance from the central spin. The only explicit dependence on k enters equation 
(13 7|) which must be averaged over using the degree distribution of nearest neighbors U(k). 
We emphasize that this approximation is exact if we expand our relations for small ?/s and 
keeping only the linear term. Note that we are finally interested in the critical behavior 
of the system where y's tend to zero and thus we expect the above approximation to work 
well in the critical region. Consequently we use the following relations to extract the critical 
behavior of Ising Model on the edge-dual of an uncorrelated random network; 



rri- 



e 2H _ e -2y 



e 2H _|_ e -2y ' 



for the magnetization, 



y, = -X)II(fc)Zn 

k 



fe-1 



J2 k ~ 1 ( ^ e {2H+y l+1 )r+K(k-l-2r-2r(k-l-r)) 



k — 1 

Z)r-o( } e {2H+y l+1 )r+k{2r-k+l-2r{k-l-r)) 

r 



for the recurrence relations defining y's and 



(42) 



(43) 



< E >= -NJJ2P(k) 



k(k-l) 



k \ 
Er=o( ) 2r ( k ~ r)e~ K2r( - k -^ +yor 
r 

Ek p—K2r(k—r)+yor 
r=0 



(44) 
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for the average of energy in the absence of magnetic field. 



At this stage it is instructive to note that we can obtain the correlation between the central 
spin and a spin at distance I by takin derivative of m with respect to Hi, the magnetic field 
acting on such a spin. To this end we need also to label the magnetic fields along with the 
y's in the above relations. Indeed in equations (}4*2|) and (|43|) the magnetic field has a similar 
index to that of y. On the other hand we have 

Ottl ~ 
Xi ■■= ^\{H l= o} = niG c {0, 1), (45) 

where hi = 2( <fc2 ^~^ :fc> ) / is the number of spins at distance I from the central spin in 
the edge-dual of random network and G c (0,l) :=< SSi > — < S >< Si >. Note that 
susceptibility is given by x — J2i Xi where from (|4*Hj) and (|4*2j) xi reads 

X ' = ^n(^) W l { ^o,. (46) 

If spins are deep enough in the network we can take all the y's equal to each other, so for 
their derivatives. After making this approximation we obtain 

dm dyi .^dyi-x i. , . 

x, = ai^ (47) 

where 



Here the index i is only to distinguish between y's in two subsequent shells and we will 
eventually set all the y's equal to each other. This quantity is determined from the fixed 
point of equation (|4*3j) . Consequently the length scale A is determined from equation 



and as expected, it will become infinite in the critical point, that is when ^ = 1- It is 
this critical behavior that gives rise to the critical behavior of % Note however that A is not 
the correlation length which is determined from the long distance behavior of G c (0, 1). 
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C. The critical behavior 



Using equation (j43}) it is not difficult to see that yo is nonzero only for temperatures less 
than T c which satisfies 

k — 1 

Er=o r ( ) e -tfc2r(k-i-r) 2 sinh(K c (2r - Jfc + 1)) 

E n W = 1- (49) 

fe — 1 

Er( ) e ^c(fc-l-2r(fc-r)) 

r 

It is a simple generalization of equation ()3ip. This equation tells us that if X c — > 0, i.e. 
at high temperatures, we have K c ~ <k'^>-<k> • ^ n °ther words the critical temperature 

P(k) is infinite. This 
231. Thus even the 



of the system becomes infinite only when the second moment o 
is the same behavior observed in uncorrelated random networks 
finite value of clustering of edge-dual networks can not significantly alter the critical point 
although T c 3> T c as we saw in the case of Bethe lattices in section IIVI 

Now let us limit ourselves to the critical region where y, H and f := | T ^ Tc I are very 

J-c 

small. We want an expansion of rh, y and < E > in terms of small deviations from the 
critical point. First note that if we change the sign of H , then by definition the sign of y 
changes too and thus the magnetization is reversed. On the other hand, energy does not 
change under this change of sign. Any expansion of these quantities in terms of H and y 
must satisfy these symmetries. We summarize these arguments in the following expansions 



m^H + y (50) 
y ^ ai {2H + y) + a z {2H + yf 
<E>^b + b 2 y 2 , 

where the coefficients a\, a 3 , 6 and 62, are found to be 



ai = l + 0{< k 4 >)f a 3 ~0(<£; 4 >) (51) 
6 ~O(<A; 2 >) 6 2 ~ 0(< A; 4 >). 
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In this equation the averages are taken with respect to the degree distribution of random 
network P(k). To be more specific let us consider a definite degree distribution, that is the 
well known scale free distribution P(k) oc &~ 7 . We consider several cases depending on the 
value of 7. 

1. The case 7 > 5 

In this case the edge-dual network behaves as a scale free network [33], with 7 = 7 — 1 > 4. 
Moreover, all the coefficients appearing in the expansions of (joTIj) are finite. It is easy to 
show using equations (|5U|l and ()51|1 that y is given by following relations 

y ~ f 3 H = (52) 

y ~ Hz f = 
y H^0,f^0. 

T 

The critical behavior of the other quantities can easily be derived from these relations 



m ~ f 2 . 



SC ~ cons., x ~ T 1 



# = 

f = 0, 



(53) 



where 5C is the change of specific heat through the critical point. Here we have only 
shown dependence of interesting quantities on f and H. Clearly these behaviors are those 
of the standard mean field model seen in the Ising model in spatial dimensions greater than 
d c = 4. This behavior is also seen in the case of Ising model on uncorrelated scale free 



random networks with 7 > 5 
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Note that due to the finiteness of all the moments of degree distribution in Bethe lattices, 
the critical behavior of Ising model on their edge-dual network also lies in this class. 



2. The case 7 = 5 



Now 7 = 4. Some of the coefficients in expansions of (|50p become infinite. To avoid 
these divergences which are an artifact of our expansion, we set a cut off for degrees which 
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is proportional to K Indeed in our expansion we used the fact that fcj/ C 1 where k is the 
degree of a nearest neighbor. Fortunately we are interested in the critical behavior where 
y — > and the the above arguments work well in that region 
Considering the above arguments, we find 



y ~ t 3 H = (54) 

. H . i 
■y ~ ( ^-)3 r = 

Tln(T) 



Thus the interesting quantities behave as 



m~r2, SC ~ in(f), x ~ — H = (55) 

rm(r) 



i 



m~( — ) 3 r = 0. 

ln(#) 



5. The case 3 < 7 < 5 



In this case the degree distribution of edge-dual network will have the exponent 2 < 7 < 4. 
Again we have to take into account the divergences appearing in the expansion coefficients. 
From equations ()5L)|) and (|5Tj) we find 



y ~ r 2 H = (56) 

y ~ H 9 t = 

T 2 



And for the thermodynamic quantities we find 



m~T2, <5C~t 7 2 5 , ^ ~ r 7 2 3 H = (57) 

~ 8-7 

m ~ 9 r = 0. 
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Magnetization 



Specific heat 



Susceptibility 



7 > 5 



cons, (cons.) 



f- 1 (r- 1 ) 



7 



T2 



y ln{H)> 



ln{f) 



Mr) 



3 < 7 < 5 



f 2 ( 7-7-3 



~ 8-7 

H— 



7-5 / S^Tl 

r 2 7-T- 3 



, _ 

t 2 r 



TABLE I: Comparison of the critical behavior of Ising model on scale free random networks 
(written inside parenthesis) and its edge-dual network. 
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We do not consider the case 7 < 3 since in this region 7 < 2 and the average number of 
neighbors is infinite in the edge-dual network although it is still finite in the corresponding 
random network. 



The above results show that these critical behaviors are very different from the ones seen in 



the uncorrelated scale free networks 



see table IJ. For example here the magnetization 



always behaves like the standard mean field case , m ~ ra, but in uncorrelated scale free 
networks this beha 
field type 0, 3- 



networks this behavior is only seen for 7 > 5 where all quantities are of the standard mean 

nn 



VI. CONCLUSION 



In summary we studied the Ising model with nearest neighbor interactions on the edge- 
dual of uncorrelated random networks. We stated a simple relation between the partition 
function of this model and that of an Ising model with next nearest neighbor interactions 
on a tree-like network. High and low temperature expansions of the partition function were 
also derived. As a simple example we studied the Ising model on the edge-dual of Bethe 
lattices using the well known recurrence relation procedure. We finally generalized this 
study to the edge-dual of uncorrelated random networks. Although the critical temperature 
of Ising model on edge-dual network is higher than the one in the random network, both 
quantities become infinite in the same point, that is when the second moment of the degree 
distribution of random network, < k 2 >, becomes infinite. This fact reflects the robustness 
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of edge-dual networks against thermal fluctuations, a property which can be attributed to 
the large number of triangles and the special structure of the edge-dual networks. We also 
derived the critical behavior of Ising model on edge-dual network of an uncorrelated random 
scale free network. The results show that this behavior is significantly different from the one 
seen in the uncorrelated random networks. 
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